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Photoassociation of spin-polarized metastable helium to the three lowest rovibrational levels of 
the J = 1, 0^ state asymptoting to 2s ^Si-|-2p^Po is studied using a second-order perturbative 
treatment of the line shifts valid for low laser intensities, and two variants of a non-perturbative 
close-coupled treatment, one based upon dressed states of the matter plus laser system, and the 
other on a modified radiative coupling which vanishes asymptotically, thus simulating experimental 
conditions. These non-perturbative treatments are valid for arbitrary laser intensities and yield 
the complete photoassociation resonance profile. Both variants give nearly identical results for the 
line shifts and widths of the resonances and show that their dependence upon laser intensity is 
very close to linear and quadratic respectively for the two lowest levels. The resonance profiles are 
superimposed upon a significant background loss, a feature for this metastable helium system not 
present in studies of photoassociation in other systems, which is due to the very shallow nature of 
the excited state 0^ potential. The results for the line shifts from the close-coupled and perturbative 
calculations agree very closely at low laser intensities. 

PACS numbers: 32.70.Jz, 34.50.Cx, 34.50.Rk, 34.20.Cf 



I. INTRODUCTION 

Photoassociation (PA) experiments in which two inter- 
acting ultracold atoms (usually ground state alkali atoms 
or metastable rare gas atoms) are resonantly excited by a 
laser to a molecular bound state provide a powerful tech- 
nique to study the dynamics of ultracold collisions P, [3l ■ 
Since the colliding atoms are so cold, the energy of the 
initial scattering state is well determined and the reso- 
nant laser energies corresponding to transitions to vari- 
ous bound states produce a very high resolution spectrum 
(< 1 MHz). At this level of precision, energy level shifts 
induced by the PA laser can be significant [1, 0, [1] and 
an understanding of the dependence of the energy level 
shifts upon the laser intensity, polarization and frequency 
is crucial. 

The short lived bound states created in photoassocia- 
tion are unique in that they can occupy both the small 
interatomic regions of conventional molecules and the 
unusual large interatomic regions of hundreds of Bohr 
radii depending upon the particular intermolecular po- 
tential. For purely long range molecules, the interaction 
between the atoms, and hence the molecular potentials, 
arises from dispersion forces that depend mainly upon 
well known atomic parameters. Extensions of the pho- 
toassociation technique can be used to determine the life- 
times of excited states 0, scattering lengths of ground 
states 01 , create ground state cold molecules [1] and drive 
other laser orientated processes. 

Although ultracold physics and the achievement of 
Bose-Einstein condensation has its roots in alkali species, 
the cooling of excited state species such as metastable 
helium opens many more opportunities for experiment. 
Metastable rare gas atoms offer exciting experimental de- 
tection strategies to study quantum gases as their large 
internal energy can be released through Penning and 



associative ionization during interatomic collisions and 
through ejection of electrons when the atoms strike a 
metal surface. Additionally, trap purity is more eas- 
ily maintained in a metastable gas, as rare gas ground 
states are not trapped by the same mechanisms as 
metastable states. The metastable atoms are gener- 
ally spin-polarized in order to suppress the autoioniza- 
tion rate and to thereby attain large numbers of trapped 
atoms. As well, no hyperfine structure is present in rare 
gas species, making them more desirable to investigate 
than many other species. 

In metastable helium a number of experimental in- 
vestigations have been conducted using photoassociation 
as the diagnostic tool. Hershbach et al. Q observed 
bound states that dissociate to the 2s'^Si-|-2p^P2 atomic 
limit, Leonard et al. [l^l studied some purely long-range 
bound states with binding energies < 1.43 GHz, disso- 
ciating to the 2s ^Si-|~2p'^Po limit, Kim et al. [ll| and 
van Rijnbach [l^ observed detailed structure of over 40 
peaks associated with bound states with binding ener- 
gies < 13.57 GHz that dissociate to the 2s^Si-|-2p3P2 
limit, and van Rijnbach [l3 | has observed six peaks ly- 
ing within 0.5 GHz of the 2s'^Si+2p^Y'i limit. Theo- 
retical analyses of the long-range bound states dissoci- 
ating to the 2s'^Si-|-2p^Po limit have been completed 
using a single channel adiabatic calculation [l^ and full 
multichannel calculations [l^. Both calculations use re- 
tarded long-range Born-Oppenheimer dispersion poten- 
tials and give excellent agreement with the measured 
binding energies. Most of the 40 peaks associated with 
the 2s'^Si-|-2p^P2 limit were identified [lB| using the ac- 
cumulated phase technique for a single channel calcula- 
tion of the bound states supported by a hybrid quintet 
potential constructed from short-range ab initio 

and ^H+,., potentials matched onto long-range retarded 



dispersion potentials. 



Recently Deguilhem et al. [if 
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have revisited the analysis of the PA peaks associated 
with the 2s '^Si+2p'^Pi 2 hmits using fully ab initio short 
range potentials. Light-induced level shifts of several 
vibrational states in the long-range J = 1,0 J poten- 
tial associated with the 2s '^Si-|-2p'^Po asymptote have 
been measured by Kim et al. [l3l and two-photon pho- 
toassociation spectroscopy has recently been used [l^ 
to accurately measure the binding energy of the least 
bound vibrational level {v = 14) of the metastable ^1]+ 
state formed during the collision of two spin-polarized 
metastable helium atoms. The measured binding energy 
Eb{v = 14) = 91.35 ± 0.06 MHz, combined with the new 
ab initio potential of Przybytek and Jeziorski p^ . 
yielded the high precision value a = 7.512 ± 0.005 nm for 
the s-wave scattering length. 

All previous theoretical investigations of photoassoci- 
ation involve limiting assumptions that may not be valid 
in the present investi gati on. Perturbative treatments of 
the radiative coupling [20l. [2l| predict a linear dependence 
of the line shifts upon laser intensity but arc only valid 
for low laser intensities. The analytical method of Si- 
moni et al. (2^ is valid for arbitrary laser intensities but 
assumes the radiative coupling vanishes asymptotically, 
thus avoiding the use of dressed states for the atoms in 
the radiation field. This is only valid when the coupling is 
negligible compared to the detuning to the atomic transi- 
tion. The most detailed treatment is that of Napolitano 
[2^ which employs a full multichannel treatment using 
dressed s- and d-wave states. However fine structure is 
neglected and it is assumed that the laser detunings are 
large compared to the radiative coupling. 

Theoretical investigations of photoassociation in ultra- 
cold metastable helium are limited to the perturbative 
analysis of the light-induced energy level shifts of the ex- 
cited J = l,0,t r ovibrational state ^ based upon the 
theories of [2ll.l23|. In this model the shifts are linear in 
the laser intensity and in the s-wave scattering length a. 
However, as the J = 1, Oj potential is very shallow, the 
laser detuning is quite small and the validity of not using 
dressed states needs investigation 

The goal of this paper is to present a complete treat- 
ment of photoassociation of spin-polarized metastable 
helium to the J ~ 1, Oj state asymptoting to 
2s^Si-h2p3po. This process is illustrated in Fig. [TJ Of 
particular interest is the dependence of the photoassocia- 
tion broadenings and level shifts upon laser intensity and 
polarization for small laser detunings comparable to the 
laser coupling. 

This paper is organised as follows. The general formal- 
ism of two atoms colliding in the presence of a laser field 
is presented in section II, a perturbative treatment of 
laser-induced energy level shifts in section III and a non- 
perturbative close-coupled treatment of the photoassoci- 
ation resonance profile in section IV. Section V presents 
our discussion and conclusions. Two appendices provide 
more details on the basis states used in our calculations 
and the evaluation of the matrix elements of the system 
Hamiltonian in this basis. Atomic units are used in the 
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FIG. 1: Photoassociation in metastable helium. Two spin- 
polarized metastable helium atoms (2s'^Si-|-2s ^Si) with very 
low kinetic energy Ek absorb a photon of frequency cj and 
form a short-lived molecule in the rovibrational level v of the 
2s '^Si -|-2p '*Po excited state. The molecule can then spon- 
taneously decay back to the disassociated metastable state (all 
other decay paths are spin- forbidden) . As illustrated here for 
the V = level, SE is the energy level shift induced by the 
laser field, AE is the laser detuning energy from the separated 
atom resonance and is the binding energy of the level. 



actual calculations, with lengths in Bohr radii oq and 
energies in Hartree Eh, = a^mgC^ = 27.211396 eV. 



II. TWO-ATOM COLLISIONS IN A LIGHT 
FIELD 

A. Hamiltonian 

The total Hamiltonian for two atoms colliding in the 
presence of a radiation field is 

H = i?mol + Hrad + -ffint (1) 

where i/moi is the total molecular Hamiltonian (in 
baryccntric coordinates) 

-f^mol = T + Hiot + Hcl + Hfs- (2) 

Here f = -{h^ /2nR){d'^/dR^)R is the kinetic energy 
operator, H^t ~ l^/i'^I^R^) is the rotational operator for 
a system with relative angular momentum I and reduced 
mass /i, Hc\ = Hi + H2 + H12 is the total electronic 
Hamiltonian of the unperturbed atoms iJi.2 and their 
electrostatic interaction H12, and His describes the fine 
structure of the atoms. The Hamiltonian for the free 
radiation field is Hrad = where a| (a^) are 

the usual creation (annihilation) operators for a photon 
of angular frequency lu^ and polarization so that the 
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field states are \n,uj^,€^) = (n!) -'^/^(a^)"|vac). The cou- 
pling between the two atoms and the radiation field is 



TO ^ — ' 



i=l,2 



where Pi = —iKVn- The vector potential is 
Mn) =Y.\-^dn) + SdnY all 



(3) 



(4) 



where £^{ri) = {h/2uj^eoVy^'^ exp(zfe • ri)e^ and V is the 
normalization volume. 



B. Close-coupled equations 

The close-coupled equations describing the interaction 
of the two atoms in an applied laser field with given an- 
gular frequency uj and given polarization are obtained 
by expanding the energy eigenstates \'^) oi H in terms 
of a basis of the form 123 



\^a,n) = \<S>a{R,q)) «) \n,uj,ex) 



(5) 



where a denotes a set of approximate quantum numbers 
describing the electronic-rotational states of the molecule 
(to be discussed in Section IIC), and q denotes the inter- 
atomic polar coordinates {9, </>) and the electronic coor- 
dinates (ri,r2). Using the expansion 



g' e' 

(6) 

where .g'(e') labels the sets of metastable (excited) states, 
in = £-1^) yields the close-coupled equations 

Y [T^,>+V,,,GAR)] +E^J"'Ge'(i?) 

= {E -nhLu)Gg{R), 



= [£;- (n- l)HGe(i?). 



Here 



T3a = -^{^a'\4^Ga{R)\^a 



(a = g or e) and 



K7 = ($e,n-l|i/i„t|$<„n) 



(7) 

(8) 
(9) 



2eoc 



{<i>e\ex-d\^„ 



where / is the laser intensity and the molecular dipole 
operator d is the sum of the atomic dipole operators 
di ~ —eri. In obtaining (fTO|) we have used [H] 



$e|-Vp,|<f3) = J(<i>e|[i?mol,d]|<I> 
TO ^-^ n 



(11) 



valid for the baryccntric frame, and have invoked the 
dipole approximation cxp(ifc • r^) w 1. This approxima- 
tion has been used in other studies of photoassociation 
SUll but does warrant some discussion. The outer turn- 
ing points of the J = 1, 0+ vibrational states considered 
here lie in the range (250-470)ao, and k = 1/3258.17 a^^ 
for the 2s ^S~2p transition so that the neglected next 
order term that would contribute is {kri^2Y ~ {kR/2)^ ^ 
5 X 10""^, comparable to other errors in our calculation. 

For photoassociation, the energy conservation relations 
arc E - nhuj = E^ + Ek and E - {n - l)huj = E^ - 
Bl + Ek + fi/\ijj where Ek = ?i^fc^/2/i is the kinetic 
energy of the colliding atoms, E"^^ are the energies of the 
asymptotically free atoms, and Alj ~ lo ~ loq is the laser 
detuning. Here huj^ ~ E^ — — E"^ is the separation 
of the unperturbed excited and ground state energies in 
the limit of zero kinetic energy. 



C. Basis states 

For two colliding atoms with orbital Li, spin Si and 
total ji angular momenta, where i = 1,2, several dif- 
ferent basis representations can be constructed. Two 
possibilities are the LS coupling scheme L = Li + L2, 
S = Si + S2 and J ~ L + S + I, and the jj coupling 
scheme ji= Li + Si, 32 = £2 + ^2, j = ji + 32 and 
J ^ j + I. The LS coupling scheme diagonalizes Hei 
whereas the jj coupling scheme diagonalizes Hfg. As the 
magnitude of the fine structure interaction is significantly 
larger than the electronic interaction at the long ranges 
of the photoassociated molecules, we use the jj coupled 
states \'yjij2jlJ'mj) (see Appendix A), where 7 repre- 
sents other relevant quantum numbers and toj labels the 
projections of the total angular momentum J on to the 
space-fixed Oz axis. A further consideration is that the 
selection rules of the laser interaction refer to the space- 
fixed reference frame and couple states of differing J and 
Toj, whereas the molecular interactions are more natu- 
rally described in the molecular reference frame in terms 
of rij, the projection of j along the inter- molecular axis 
OZ. Hence we choose the hybrid Hund case (c) molecular 
basis 



\<^a{R,q)) = \lll2jlj2j^jW,Jm,j) 



(12) 



(10) 



where w is the symmetry under inversion of the elec- 
tronic wavefunction through the centre of charge (see 
Appendix A). 

The matrix elements of the various contributions to 
the Hamiltonian in this basis are derived in Appendix B. 
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III. PERTURBATIVE TREATMENT 
A. F-operator technique 

For low laser intensities Hint can be treated as a 
perturbation. The unperturbed Hamiltonian is then 
Ho = i?moi + ^rad and has eigenstates of the form 
R~^G^{R)\ip^) (g) \n,uj, ex) where the unperturbed radial 
functions G^(i?) satisfy Q with V^f omitted and E re- 
placed by the total unperturbed energy Eq. 

To calculate the second-order laser-induced energy 

(2) 

level shifts 6Ee of the excited bound states |e) we note 
that [l^ have shown that the adiabatic approximation 
has little effect on the bound levels, making it possible to 
consider each excited state as an isolated state. The adi- 
abatic states are formed by diagonalizing iJmoi without 
the radial kinetic term and can be written as 



CeaiR)\^a) 



(13) 



to eval- 



where Ealt^ea(i?)P = l. 

We employ the F-operator technique [23, [2? 

(2) 

uate SEe ■ This approach has recently been used by 
Beams et al. [2^, to treat spin-dipolc interactions as 
a perturbation and we adapt their formalism here to the 
laser- matter interaction Hint • The shift is given by 



(e|i/i„tF|e) 
(elHintFle) 



(e|F|e)(e|i/int|e) 



(14) 



where the second term is zero as (e|ifint|e) vanishes due 
to dipolc selection rules. The operator F satisfies 



[F,Ho]|e) =Hi, 



(15) 



We assume the states , n)} and 1-0°, n—1) = iV'e)® 
\n — l^Lo^ex) form a basis and expand F\e) in terms of 
them: 



Fie) = R- 



a' 



so that p5)) gives the coupled equations 



2fj.dR^ 
and 



[So - (n - 1)M - K^'^ - 0,|(i?) 



(16) 

feiR) = 

(17) 



{Eq — nTiLo) 



fAR) 



_2fi dR^ 

E + K's] MR) = G",{R)V^J^HR)- (18) 



We have assumed no i?— dependence of the basis states 
|$a) but include the correction term V?^*-^ for the excited 



state (see Appendix B). The matrix elements are 

= ($,,|7?,ot|<I>3) (19) 

and 0+(i?) is the adiabatic potential of the excited state: 

OtiR) = (V'°l[4ot + ^^ei + HtM)- (20) 
The interaction term is 

Fj"*(i?) = (<i>„n|A„t|V',",n-l) (21) 

and G'g(i?) is the radial eigenfunction for the rovibra- 
tional level v of the excited state. 

Substitution of the expansion (fTB)l into (fH|) yields 

^^i'^ = E /^^ G"ARK''\R)m)- (22) 

Consequently we need only solve the inhomogeneous 
equations ([T5)) . These are solved using a DVR method 
[31| with a cosine Fourier basis. As the laser interaction 
region requires closely spaced grid points and the asymp- 
totic region very few points, a scaled grid R = ({t) [20| 
was used where the mesh in t is equispaccd. A quartic 
scaling was chosen as it does not modify the boundary 
conditions. 

B. Application to metastable helium 

We consider photoassociation to an excited long-range 
bound level in the J = 1, mj = 1, adiabatic potential 
that asymptotes to a He(2s '^Si) + He(2p'^Po) diatom 
(see Fig. [1]). The adiabatic state can be expressed as a 
combination of the four basis states | j2 , j) 



|0,1> 



11,1) 



|2,1) |2,3) 



(23) 



where we have suppressed the common quantum numbers 
71 ^ (Li,5i) = (0,1), 72 = (7.2,52) = (1,1), ji = 1, 



0, J 



771 J 



1 and u) = 1. Coriolis couplings 



and the off-diagonal kinetic terms (sec Appendix B) have 
been ignored, as their effect has been shown to be small 
at long range [T3 |. 

An analysis of iJjnt shows that only a small mani- 
fold of metastable basis states are coupled to the excited 
state. Using the reduced state notation of J77ij), 
the states coupled by cr" laser light are 

100,22), |20,22), |21,22), |2-1,22), 

(|22,22), |2-2,22)) (24) 

whereas the states coupled by (7+ laser light are 

|00,00), |00,20), |20,00), |20,20), 
|21,10), |2-1,10), |21,20), |2-1,20), 

(|20,10), |22,20), |2-2,20)). (25) 
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TABLE I: Perturbative line shifts per laser intensity in units 
of MHz/(W.cm^^) for the long range bound states of the 
He(2s"^S'i) + He(2p''Pi), 0^ configuration. Results are given 
for calculations without and with the correction V?^. 



Level 


Polarization 


No V^"" 


With <^ 


Ref [5] 


w = 




-6.439 


-6.507 


-6.37 




a+ 


-7.724 


-7.784 


-7.36 


V = 1 




-11.662 


-11.748 


-11.70 




a+ 


-10.205 


-10.270 


-10.25 


v = 2 




-29.442 


-29.692 


-29.57 




a+ 


-24.675 


-24.877 


-24.11 



The bracketed states are coupled by the Coriolis terms 
of iJrot to other metastable states and hence are 
indirectly coupled to the excited state. The im- 
plicit quantum numbers for these metastable states are 
71 = (Li, ^i) = (0, 1), 72 = (i2, ^2) = (0, 1), ii ^32^1 
and w = 0. 

The calculations require Born-Oppenheimer potentials 
for the and states of the two metastable 

He(2s'^Si) atoms and the J = IjOj potential of the 
excited He(2s^Si) -I- He(2p^Po) system. For the quin- 
tet potential we use that given by Przybytek and 
Jeziorski which includes adiabatic, relativistic and 
QED corrections. For the singlet potential we use 
[23] a potential constructed from the short-range Miillcr 
et al [32] potential exponentially damped onto the quin- 
tet potential at long range. The excited state poten- 
tials are constructed from the twelve Born-Oppenheimer 
dispersion potentials [H fsAiR/XjCsA/R^ + Cqa/R^ + 



C^^/R^, where /3A is an R- and A-dependent retardation 
correction [s^, A = 3258. 17ao, where A = A/(27r) and A 
is the 2s '^S — 2p transition wavelength. In particular, 
the excited state |0+) is a linear combination of the ^5]+, 
^n„, ^Uu, and states. The C3S coefficient is ±20^ 
and Can is iCa, where C3 = 6.41022 EhUQ and is propor- 
tional to the square of the 2s^S — 2p^P transition dipole 
matrix clement. For the van dcr Waals coefficients we use 
Ces = 2620.76 Ehal, and Cm = 18^6.60 EhU^. The C^^ 
terms are [H C+v = 151383 Ena^ C^^ = 297215.9 i;;jag, 



97244.75 S^ag and C, 



7gn = 162763.8 £:ftag. Here 
the superscript indicates the sign of (—1)'^+"' where S is 
the spin of the state and w ~ 0(1) for gerade (ungerade) 
symmetry. The C^/R^ term is the dominant contribu- 



tion to the dispersion potential for the purely-long-range 
states. Also required is the 2s —^ 2p atomic dipole mo- 
ment = 2.14583 X 10-29 c.m and the fine structure 
sphttings. For He the 2p^Po and 2p^Pi levels lie 31.9081 
GHz and 2.2912 GHz above the 2p^P2 level, respectively. 

The results for the three lowest vibrational bound lev- 
els and for laser polarizations cr^ are given in tablep) In 
order to compare with the results of Portier et al. Q, we 
have used Ao; = 0, corresponding to a zero kinetic en- 
ergy. There is good general agreement between the two 
sets of results with the small differences probably due to 
the potentials used. 

IV. CLOSE-COUPLED CALCULATION OF 
PHOTOASSOCIATION PROFILE 

A. Hamiltonian matrix 

The close-coupled calculation is non-perturbative in 
that the differential equations ([7]) are solved without as- 
suming the laser interaction is weak. In this approach the 
scattering matrix elements are calculated and the PA pro- 
file obtained for various laser detunings, intensities and 
temperatures. 

Although the method presented here is quite general, 
for ease of visualization we explicitly formulate it for pho- 
toassociation in metastable helium of the subset of states 
coupled by a~ polarized light in (|24p . Specifically these 
states la) arc 



1) 


- \J 


= 0, nj = 0, J = 


2, mj = 


2) «) \n 






2) 


= 1.? 


= 2, % = -2, J 


= 2, mj 


= 2) (g) 


n, 




3) 


- b- 


= 2,% = -1, J 


= 2, mj 


= 2) «) 


n. 




4) 


- 1.? 


= 2, % = 0, J 


= 2, mj 


= 2) (g) 


n. 




5) 


= \3 


= 2,%= i,j 


= 2,m.7 


= 2) ® 


n, 




6) 


= \j 


= 2, % = 2, J 


= 2,m.7 


= 2) ® 


n, 




7) 




^J=l,m,, = l 


)(E)\n- 


1,^,£\) 







(26) 



for which the close-coupled equations ([7]) reduce to 



E 



2//di?2 

where the potential matrix W is 



G„'(i?) = (27) 



W = 



- Ek 






















^5"' - Ek 





















y,'' Ek 


L^i 


















~ Ek 


£01 
















LOI 


- Ek 


£12 
















£12 


Vi - Ek 












nnl 




K- 





(28) 
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Here 



where the dresscd-state radial amphtudcs satisfy 



v^%, = + {jn,\H,,,\jn,), (29) 

is the effective ground state potential, comprising the 
Born-Oppcnheimer potential for the state and 

the diagonal rotational couplings, Ve = 0+(i?) — iT/2 is 
the adiabatic potential of the excited state, and AE = 
TiAlu + Ek — B1 . The off-diagonal terms are the Coriolis 
couplings 



i''-''^ =(.? = 2,^-|i?,ot|j = 2,r!, 



(30) 



V'"'(i?) expressed 



and the radiative couplings 

in terms of the Rabi frequencies {^'^j^^^. 

The decay of the excited state through spontaneous 
emission is represented by the molecular width V = 
r(i?/A, Fat) where the i?-dependence arises from the re- 
tarded dipole interaction expressed in the Hund's case (a) 
states |33j and its subsequent diagonalization as part of 
the formation of the adiabatic potential. For our unger- 
ade system, F « 2Fat for the interaction region and 
F K, Fat for the asymptotic region R ^ X. The spon- 
taneous emission atomic line width for the triplet he- 
lium 2s-2p transition is Fat = 1-626 MHz. Note that the 
asymptotic energies E^^ cancel in (pS)) . 

The definition of scattering matrix elements requires 
asymptotically free states. However the asymptotic form 
Woo of W is not diagonal as the Rabi couplings do not 
vanish at large R and the separated atoms remained cou- 
pled by the laser interaction. Two options are available 
to determine the S'-matrix elements: (i) transformation 
from the basis states \a) to the dressed-state basis states 
El 1/3) = Ea Uf3a.\a) in which WD = U-'W^U is di- 
agonal or, (ii) introduction of a modified radiative cou- 
pling for R greater than some large value Rz'- 



v:f{R) 



for R< R. 



V^\R) exp [-p{R - i?^)2], for R> R^ 



(31) 

This coupling vanishes asymptotically and simulates the 
experimental conditions of the laser being switched off 
before and after the experiment. 



B. Dressed States 

Napolitano [2^ has considered the case of six dressed 
states with laser couplings restricted to two pairs of states 
and has obtained analytical results valid for large red 
detunings AE ^ F and AE ^ fig ^ . As we wish to avoid 
such assumptions and also develop a procedure valid for a 
larger number of coupled states, we use direct numerical 
diagonalization. 

In terms of dressed states the expansion of the state 
vector becomes 

= E Gc.{R)\a) = R-' Gp{Rm (32) 



E 



h' dp 

2^ rfi?2 



Gp,{R)^Q (33) 



and = U ^WU. Asymptotically the dressed states 
decouple and satisfy 



dR^ 



2^lEp 



Gp{R) = 



(34) 



where the energies i?/? = — W^(i? = oo) can be complex. 
Defining fc^ = yJ^p.Epj'h = k'^ + ik^, the asymptotic 
solutions of ((34|) have the form 



GpiR) ~ cie 



(35) 



The open-channel and closed-channel manifolds are then 
identified respectively as those states that persist or van- 
ish as i? — > oo. Therefore we require fc^ to be zero for 
the open channels. 

After diagonalisation we find that two of the energies 
Ep are nondegenerate and their associated dressed states 
include undressed excited state contributions. The third 
value of Ep is purely real, five-fold degenerate and mixes 
only the undressed open channels for the ground states. 
As the diagonalization is numerical, the particular com- 
bination of ground states in this degenerate subspace is 
dependant upon the numerical procedure used and the 
combinations will not vary smoothly as the detuning is 
changed. However, the S'-matrix elements that involve 
only the states with non-degenerate energies will be guar- 
anteed to vary smoothly. 

For the orthogonal polarization we briefiy note that 
the above analysis applied to the 11 metastable states 
listed in ([25|l exhibits very similar behaviour. Three 
unique values for Ep are observed, two that are non- 
degenerate and one that is 10-fold degenerate. No differ- 
ences in the analysis or numerical implementation apply 
other than to include a larger number of states in the 
system. 



C. Modified Radiation Coupling 

Because the dressed states introduce complex asymp- 
totic energies, the simplicity offered by using a modified 
radiative coupling becomes more attractive. To imple- 
ment the modified coupling, we use the undressed form 
of the state vector (1321 to obtain 



E 



d^ 

' 2ii dR^ 



Ga'(i?) = (36) 



where W^^, is identical to Waa' except that V^g^ is re- 
placed by V"g^. This system consists of six open chan- 
nels, all with identical values of ka 
closed channel. 



\/2fiEK /h, and one 
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D. Extraction of S'-matrix elements 

The S'-matrix is determined by matching the asymp- 
totic solutions of ([551) o^' (EH) to the combination [S^l 



G = H^ A + HlB = (H^_ - HlS)A 



(37) 



where G^~fi is the matrix of solutions formed from G-y{R) 
with the second subscript 7' labelling the linearly in- 
dependent solutions generated by different choices of 
boundary conditions. We have introduced the notation 
7 = /3 for the dressed states approach and j = a for 
the modified coupling approach. The diagonal matrices 



Sjy'h^ for the open channel scattering states 



have the asymptotic form of outward and inward travel- 
ling waves: 



R — >OD 



(2|M)- 



(38) 



whereas for the closed channels the asymptotic forms are 



R — >oo 



(2 m 



(39) 



The matrices A, B, and S = BA ^ have the struc- 
ture 



4 4 

^^00 ^oc 

4 4 



(40) 



where the labels o and c refer respectively to open and 
closed channels. We require only the open-open contri- 
butions Soo which, as we can enforce Aco = by care- 
ful choice of boundary conditions, are given by Soo = 
— R 4 1 

As Ej3 ^ Ep + iEJ^ and fc^ can be complex for the 
dressed states, the asymptotic solutions do not have the 
usual plane wave oscillatory form and care must be taken 
during the matching process. In the present case Ep is 
complex for only two of the dressed states. One is clearly 
a closed channel as it has Ep < 0. The other has EJ^ > 
and a relatively small value of i?^. Without the imagi- 
nary component, this channel would have a purely real 
value of kfj and be considered an open channel. However, 
as kf3 is complex, the state has the asymptotic form psp , 
consisting of exponentially increasing and decreasing so- 
lutions. To treat the solution rigorously, we must enforce 
finiteness by discarding the exponentially increasing so- 
lution and integrate out to a sufficiently large value of R 
that the exponentially decreasing solution has completely 
dampened, indicating that the channel is closed. The 
S'-matrix can then be created from the remaining open 
channels that have purely real Efj components. The dif- 
ficulty with such an approach is that ^ fc^, requiring 
the integration range to be very large, imposing a heavy 
computational burden. Note that kp and k^ will also 
vary with the laser detuning and laser intensity, further 
complicating the problem. 

An alternative approach is to consider the problematic 
channel as a pseudo-open channel. By redefining the S- 
matrix, such that the asymptotic functions are matched 



at a finite distance R = i?max instead of i? — > c», we can 
ensure the pseudo-open channel has a finite inward and 
outward flux. This is performed by treating the terms 
gTk},R approximately constant in the local region of 
-Rmax and matching to the oscillatory behavior of e^''^^'^. 
This type of matching obviously has a dependence upon 
the matching point but we expect it to not vary the shape 
of the profile significantly if i?max is chosen outside the 
interaction region. The resultant S'-matrix element Sfjip 
then gives the probability that the system with unit fiux 
in an incoming channel |/3) at i?max makes a transition 
to an outgoing channel |/3') at Rmax- We note that this 
choice of i?max bears some similarity to the choice of Rz 
for the modified coupling method. 



E. Detection of the Resonance 

The photoassociation resonance can be studied by 
analysing the loss from the excited state due to sponta- 
neous emission so that the photoassociation lineshape is 
due to the emitted photons and therefore proportional 
to the loss of unitarity of the S-matrix. For atoms 
colliding in the entrance channel I7). the loss rate is 



-7 — \ 
sion is 



11-^0-^'^°'°") where the cross section for photon cmis- 



photon 
7 



1,2 
'"7 



E 

7' 



(41) 



and (• • • ) denotes a thermal average over a distribution 
of the asymptotic relative velocities Vj = hky/jj of the 
two colliding atoms. For temperatures of order 1 /iK this 
thermal averaging can be ignored. The energy depen- 
dence of crP'^°*°" is that of a peak superimposed upon a 
slowly varying background and can be well fitted in the 
region of the peak by a Fano profile HI] of the form 



^photon ^ _ 



' l-fe2 



(42) 



where Abg(e) describes a linear background, A-^^s is a con- 
stant, e = (-E — ii'rcs)/(rres/2) IS a normalized energy, and 
-Eros and Fres are the position and full width of the res- 
onance. The Fano parameter g is a measure of the ratio 
of the direct (background) to resonant scattering. As we 
often deal with Lorcntzian-like behaviour that occurs in 
the limit q 00 and A^^s ^ 0, it is simpler numerically 
to instead match to the form 



(7. 



photon 
7 



(43) 



where p = l/q and A[.^g ~ (Z^rcs- For the present cal- 
culations we find p < 10~^ except for the dressed state 
profiles at high intensity where p 0.5. 
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F. Numerical issues 

To solve the dressed problem, the undressed asymp- 
totic potential matrix Woo is diagonalized by a Her- 
mitian eigendecomposition. The renormalized Numerov 
[stI boundary value method is then used to integrate ei- 
ther the dressed equations ((33|) or the modified coupling 
equations ([36|) . The outer boundary used is either i?rnax 
for the dressed states or, for the modified coupling, a 
value greater than such that the radiation coupling 
has been completely turned ofl^. To ensure the asymp- 
totic solutions are uncoupled, we choose i?niax and to 
be greater than 10'* oq. Typically we use Rz = 2x lO^ao 
and p ~ 0{5k~^) = 4 x 10~^aQ ^. The linearly indepen- 
dent solutions Gj^', 7' = 1, — 1, are generated by 
choosing A/"— 1 linearly independent boundary conditions. 

A kinetic energy of 10^** Eh (2.1 fiK) is chosen, which 
places the system just above temperatures for which re- 
cent experiments have reported Bose-Einstein condensa- 
tion and which docs not introduce a prohibitively large 
outer boundary for the numerical integration. Note that 
quantitatively, at an intensity of 0.7 W/cm^, the matrix 
elements in W^o are O(10~^*) and O(10~^) for the diag- 
onal elements Ek and AE respectively and O(10~^) for 
the off diagonal Rabi couplings f^2/+i- 



G. Results 



20|-^ 
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FIG. 2: PA cross section profile calculated using dressed 
states. Results shown are for polarization and a low in- 
tensity of 64 mW/cm^. 
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We first consider the results for low to moderate laser 
intensities as they exhibit behavior similar to regular 
spectroscopic profiles. As the laser intensity is increased, 
unusual aspects of the dressed and modified coupling be- 
come apparent and we shall discuss these separately. 

The central positions and broadenings of the dressed 
PA profiles were determined from the fits to the Fano 
profiles 142]) for the cross section fjP'^°'°" as a function 
of laser energy. For nondegenerate channels such as 
the pseudo-open channel in the present investigation, the 
variation of (t^'^°*°" with laser energy is smooth. How- 
ever, as previously mentioned, the numerical diagonal- 
ization process arbitrarily selects the degenerate dressed 
states {|/3(i)} in the subspacc £d that they span. This 
means that crj^^"*"" for any (3 ^ (3d will not vary smoothly 
with laser energy. 

Fortunately, for these degenerate dressed channels we 
are permitted to analyse various combinations of (t^'^°*°" 
as the physical behaviour of the system cannot depend 
on the choice of basis in Sd- The simplest choice is an 
average of all channels, i.e. 



photon 
Mresscd 



1 



rid 



^ ^ ^photon 

13 



(44) 



where nd is the number of degenerate states. Results 
for the case of a~ polarization and a low intensity of 64 
mW/cm^ arc shown in Fig. [21 
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FIG. 3: PA cross section profile calculated using the modified 
coupling. Results shown are for a~ polarization and a low 
intensity of 64 mW/cm^. 



For the modified coupling method, the cross section 
^photon -j-Qj. ^j^j, undressed channel \a) is more straight 
forward to analyse in terms of experimental conditions. 
We choose to form a similar quantity 



photon 



1 



'^modified ^ 



photon 



(45) 



where rio is the number of open channels, and again fit 
the resonances to Fano profiles. In this way we can com- 
pare the profiles with those of the dressed state calcula- 
tion. Results for a~ polarization and an intensity of 64 
mW/cm^ arc shown in Fig. [31 

The two profiles shown in Figs [2l and [3l are very similar 
apart from the suppression of the background at small 
dctunings for the dressed state profile. Importantly, the 
resonance parameters Ej-es and Fres obtained from fits to 
the two profiles are identical. The background loss be- 
havior present in the profiles is due to the dominance of 
the off-diagonal Rabi couplings over the diagonal terms 
AE and Ek of the potential matrix for the asymptotic re- 
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FIG. 4: PA cross section profile calculated using dressed 
states. Results shown are for polarization and a high 
intensity of 2.6 W/cm^. 



FIG. 6: (Color online) Dependence upon laser intensity of the 
line shift of the v = 0,1,2 levels for polarization coupled 
to the 0^ state. 
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FIG. 5: PA cross section profile calculated using the modified 
coupling. Results shown are for polarization and a high 
intensity of 2.6 W/cm^. 



FIG. 7: (Color online) Dependence upon laser intensity of the 
line width of the u = 0, 1, 2 levels for polarization coupled 
to the Ou state. 



gion. This situation arises because of the shallow nature 
of the excited state potential and the ultracold temper- 
ature used. If the well is artificially deepened, or much 
larger kinetic energies are used, then this background 
completely disappears. 

The PA profiles for high laser intensity show some un- 
usual behavior. Results for the dressed state and mod- 
ified coupling profiles at an intensity of 2.6 W/cm^ are 
shown in Figs [3] and [5] respectively. These spectra ex- 
hibit two features not present in previously calculated PA 
profiles; strong interference between the resonance and 
background contributions and a reduction in the overall 
magnitude of the dressed profile. The interference fea- 
ture is apparent in the magnitude of the corresponding 
Fano parameter p = !/<? ~ 0.5, indicating that these res- 
onances exhibit a severe departure from Lorcntzian-like 
behavior. 

The dependence of the resonance position and width 
upon laser intensity is shown in Figs [5] and [7] respec- 
tively for the case of cr" polarization. Similar behavior 



is obtained for (7+ polarization. The intensity depen- 
dence of the line shift and width is very close to linear 
and quadratic respectively for v = 0,1 but departures 
from these dependencies are evident for v = 2. This can 
be seen from the fit parameters given in Table [TT] which 
shows that the coefficients S2 of the quadratic correction 
and W3 of the cubic correction to the intensity depen- 
dence of the shift and width respectively for v = 0,1 are 
very small. 



V. DISCUSSION AND CONCLUSIONS 

We have investigated the line shifts and widths for 
photoassociation of spin-polarized metastablc helium to 
the three lowest rovibrational levels of the J = 1, 
state asymptoting to 2s ■^Si-|-2p^Po using two variants of 
a non-perturbative close-coupled treatment, one based 
upon dressed states and the other on a modified radia- 
tive coupling which vanishes asymptotically. We have 
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TABLE II: Parameters for a quadratic fit sil + 82!^ and cubic fit too + wil + wzl'^ + ^3/'' to the dependence of the line 
shifts and the line full widths respectively (in MHz) upon laser intensity / (W/cm^). The value /max denotes the approximate 



maximum intensity up to which the resonance 


peak is 


clearly discernable and there 


is no overlap with 


neighboring peaks 




Level 


Polarization 


Shift 








Width 




-^max 






Sl 


S2 


Wo 




W2 


UI3 




w = 




-6.514 


-0.023 


3.233 


0.441 


0.304 


-0.0093 


7.0 


w = 


a+ 


-7.781 


-0.021 


3.226 


0.730 


0.329 


-0.0091 


7.0 


v = l 


a~ 


-11.77 


0.074 


3.231 


1.51 


2.87 


-0.34 


3.2 


v = l 


cr+ 


-10.30 


-0.018 


3.219 


1.47 


1.70 


-0.10 


3.2 


v^2 


a~ 


-29.79 


3.90 


3.216 


6.25 


36.9 


-21.9 


0.4 


v = 2 


a+ 


-24.95 


1.13 


3.197 


5.44 


18.2 


-4.20 


0.6 



also calculated the shifts using a second order perturba- 
tive treatment. 

The main physical interest is in the properties of the 
PA resonance profiles and our results for the shifts and 
widths of these profiles obtained using the two variants 
of the non-perturbative calculation agree to better than 
2%. Both methods indicate that there is significant back- 
ground loss for this metastable helium system, a feature 
not present in studies of PA in other systems, which is 
due to the shallow nature of the excited state potential. 

The behavior of the PA profiles for high laser inten- 
sity is quite different for the dressed state and modified 
coupling methods. Physically, the behavior of the mod- 
ified coupling profile appears to make sense; if the laser 
is intense enough that most loss occurs outside the colli- 
sion region, then the particular resonance of the excited 
state will have little effect on the profile. Hence a form 
of saturation is observed. Note that this first occurs in 
the higher vibrational levels as the intensity is increased. 
The behavior of the dressed states profile is very different 
in that the overall magnitude of the profile decreases at 
higher intensities and lower detunings. It is difficult to 
relate this behavior to the collision process without the 
appropriate association of the dressed state description 
with the experimental, undressed states. When the laser 
is switched on, the initial undressed states must be trans- 
formed to the dressed states and in this process the un- 
dressed open channels acquire significant components of 
the pseudo-open and closed dressed channels and thereby 
suffer loss. In the limit of infinite intensity, the open 
dressed channels become completely uncoupled from the 
pseudo-open and closed channels, and the loss from the 
dressed states is entirely due to the activation of the laser. 
This suggests that the general interpretation of quanti- 
ties constructed from dressed S-matrix elements at large 
intensities requires a proper description of the activation 
of the laser in the formalism. Fortunately, this detail is 
not required to obtain the resonance parameters as both 
treatments result in nearly identical shift and width val- 
ues. 

The results for the line shifts from the close-coupled 
and perturbative calculations agree very closely at low 
laser intensities. The si values in Table |TT] differ by 
^ 0.5% from our perturbative results (see Table |T| cal- 



culated with the correction Vg included. The small 
differences can be explained by the fact that a finite ki- 
netic energy (2.1 fiK) was used in the close-coupled cal- 
culation whereas a zero kinetic energy was used in the 
perturbative calculation. We note that an increase in 
the kinetic energy by an order of magnitude decreases 
the non-perturbative line shift fit parameter by approxi- 
mately 2%. 

At higher intensities the perturbative results remain a 
fair approximation to the non-perturbative results, how- 
ever small non-linear differences are evident, affecting the 
line shifts by up to 5%. Additionally, a maximum inten- 
sity is found that limits the visibility of the resonances. 
This laser saturation should be observable in the labo- 
ratory and suggests that ultracold photoassociation to 
shallow potentials is only useful at lower laser intensi- 
ties. Although not reported here, it can be shown Q 
that the perturbative line width is also linearly depen- 
dent upon laser intensity. This is obviously not true in 
the present non-perturbative calculation, where a signif- 
icant quadratic behavior is observed. 

We conclude by noting that our theoretical findings are 
consistent with the limited experimental data available 
[l3|. Our values for the shift ratios si(w = l)/si{v = 
0) and si{v = 2)/si{v = 0) of 1.81 and 4.57 for ct" 
polarization lie within the experimental values of 1.71 ± 
0.14 and 4.20±0.48 respectively, and our v = widths of 
3.24 MHz at 9 mW/cm^ and 11.9 MHz at 5 W/cm^ are 
comparable to the experimental values estimated from 
Fig. 2 of [l3| of 3 MHz and 10 MHz respectively. 
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APPENDIX A: THE \j1j2jnj, Jmjp) BASIS 

In the jj coupling scheme, the total angular momenta 
ji of each atom are coupled to form the total electronic 
angular momentum j = ji -\- J2 which is then coupled to 
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the relative angular momentum I of the nuclei to form 
the total angular momentum J = j + I. This gives the 
basis states in the space-fixed reference frame 



E 

nijini 



Cr "m,mj \ljlhjmj)\lmi) 



(Al) 



a plane containing the intermolecular axis. For identical 
nuclei the symmetric states for jj and LS couplings are 

+ (-l)^-|72j27iJijf^j» (A7) 



and 



where \jj1j2jTnj) SlTC the eigenstates of the two asymp- 
totically free atoms and \lmi) = Yi^rniiS,4') are the rel- 
ative motion eigenstates. Here 7 = (71,72) and 7^ = 
{7;, Li, Si} where Li and Si are the total orbital and to- 
tal spin angular momentum respectively of the individ- 
ual atoms. The label 7^ denotes any additional quantum 
numbers needed, including those specifying the electron 
configurations. All subscripted m quantities denote pro- 
jections on the Oz axis of the space-fixed frame. 

The transformation from the space-fixed frame to the 
molecular frame is 



= NLs{\iii2LsnLns) 

+ {-ir^'h2iiLsnLns)) (A8) 



where the explicit ordering of 71 and 72 indicates the 
order in which the angular momenta are coupled and 
p, which labels the symmetry under permutation of the 
labels 1 ^ 2, is uniquely related to w. For the jj state 
(|A7P , pjj — w + wi + W2 + N + ji + j2 — j where Wi is the 
symmetry under inversion of the electronic wavefunction 
of atom i about its nucleus and N is the total number 



\lJ1J2jmj) 



D: 



of electrons per atom. The LS state (|A8p has pls ^ 
*^ j^^. ((/), e, 0)\-fjij2jnj) (A2) w + wi + W2 + N + Si + S2 - S + Li + L2 ~ L. 

The normalization constants Njj and N^s are 



where subscripted Q quantities indicate projections along 
the intermolecular axis and D^^ ,^, is the Wigner rotation 
matrix (ssj . Expressing the relative motion state \lmi) 
as a rotation matrix gives 



\lj1j2jlJmj) = ^ Ci 



u 

mjmiTiij 



V 



21 + 1 



4tt 



-i?:;o(^,e,o) 



X Di:^^u,0,o)\jnj2jn,). 



(A3) 



Combining the rotation matrices and using the sum rule 
for Clebscli-Gordan coefficients reduces (IA3I) to 



\lJ1J2jlJmj) = ^(-l)^"-^^^C^f_o,o 

X Ni^^^{9,cl,)\jnj2jn,) (A4) 
where .0 . is the symmetric top function defined as 



rn J £1 j 



2J+\ 

47r 



(A5) 



Equation (|A4[) can be interpreted as a coupling of j and 
J to result in / and naturally introduces the basis states 



N,. 



75 for (71, ji) ^ (72,^2) 
5 for (71, ji) = (72, j2) 



and 



LS 



for 71 ^ 72 
2 for 71 



72 



The transformation between the two bases 
(IA8I1 is 



(A9) 

(AlO) 
3 and 



|7lil72j2.)^^jW) 



E 



N. 



33 Tjhhj^j 



LS 



LSni^Us 
X \ji^2LSQl^s'w) 



where 



LSUl^s 



[(25+l)(2L+l)(2ji + l)(2j2 




(All) 

l)]i 
(A12) 



where the {• • •} is the Wigner 9— j symbol and the implicit 
set of quantum numbers (71, 72) has been suppressed. 



(A6) 



APPENDIX B: MATRIX ELEMENTS 



The derivation of (|A6|) has not taken into account any 
of the symmetry requirements of the system. Following 
[3^, we define symmetric states that are eigenstates of 
the operator /, the inversion operator of the total wave- 
function through the centre of charge of the molecule. 
The eigenvalues of / are (— 1)*" with u; = for gerade 
symmetry and w ~ 1 for ungerade symmetry. For rij = 
states, we can also identify the quantum number of (T„, 
the refiection operator of the total wavefunction through 



1. Kinetic terms 



The radial kinetic term has the form 



2^iR V di?2 '^^ dR^^UR 



Ga{a' 



dR^ 



(Bl) 
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where \a) = |$a(-R^<?)) and a represents the quan- 
tum numbers {71, 72, Ji, J2j J, ^j, w, J, mj}. As the basis 
states are assumed to vary Uttle with respect to i?, the 
last two terms in (|Bip are neghgible at the long ranges 
considered in this investigation. 

In this investigation the adiabatic excited state ji/'g) 
is a combination of these basis states (see p3)) ) and the 
coefBcients Cea{R) do vary considerably with R. If we 
assume no other excited states are coupled to the system 
(i.e. the adiabatic approximation is exact) then we only 
require the radial term (jBip with \a) = \a') = The 
second term of (IBll) is then zero as 



Electronic term 



The third term, however, is non-zero and gives rise to the 
kinetic correction term 



V.^'^=^-S^j:C.aiR)^^^%^. (B3) 



2fiR 



dR^ 



The matrix elements of the rotational kinetic term 



2/ii?2 



(B4) 



are evaluated using the expansion of r in terms of ladder 
operators 



l'' = J^+ ? - (2 J.i. + J+]- + J-J+) 



(B5) 



where the subscripts refer to molecule-fixed axes, J± = 
Jx ± iJy, and j± = jx ± ijy The action of J± is irregular 
[iot due to the rotation of N^^^^, and is given by 



Hence the matrix elements of r are 

{a'\P\a) = hHpp, { [J( J + 1) + j(j + 1) - 2^2] 5n'p, 

- Kjjn, ^ ,0,-1 - Kjjn, ^ ,n, +1 } (B7) 

where \a) = |<f>a(i?, g)), a represents the quantum num- 
bers {71, 72, Ji,j2,j, f^j, J, mj}, p denotes the set of 
quantum numbers {71 , 72 , Ji , ^2 , i, w , J, mj } and 

= [J(J + l)-r!,(%±l)]^ 

X [.7(j + l)-f^,(^±l)]^ (B8) 

The terms non-diagonal in Clj are called the Coriolis cou- 
plings and are often negligible. This is the case for purely 
long-range bound states in metastablc helium (l^ . 



We wish to express the matrix elements of iJ^i in terms 
the Born-Oppenheimer potentials '^^^^A^{R) defined by 
the eigenvalue equation 



H^ilLSnLftsw) 



r2S+i A± 



A^iR} + E°°]\LSnLnsw) (B9) 



where A = I^IlI and E°° is the asymptotic energy of the 
state. 

The matrix elements in the basis (|A7[) are evaluated 
by transforming to the basis (|A8|) using (|Alip and then 
applying (|B9[) to obtain 



(a' I Hell a) 



E 



m2 , , , 



■^m' ^ LsnL^s 

p«+iA±(i?)+£r]^i\^a 



(BIO) 



where 77 = {71, 72, fij, w, J, mj}. For our case of 
metastablc helium involving the 2s2s and 2s2p states. 



3. Fine-structure term 

The total fine structure term Hfs is the sum of the fine 
structure terms for the individual atoms and is diagonal 
in the basis (jA6p in the asymptotic region of free atoms. 
For the long-range molecular states considered in this 
work, the total atomic angular momenta are considered 
to be approximately good quantum numbers. Hence the 
matrix elements of the total fine-structure term are 



{a'\Hfs\a) 



Saa'iAE. 



fs 

71 Jl 



where iJj'g represents the fine-structure interaction and 
AE^.j. the fine-structure splittings for atom i. 



4. Laser interaction term 

For radiation of a given circular polarization €\ in the 
space-fixed frame, where A = 0, ±1 for tt, ct* polarization, 
the laser-matter interaction (jlOp can be expanded in a 
spherical basis using 



E 

C=o,±i 



(B12) 



where (ex)-^ = (5a, 5. The matrix elements of ffmt are, 
after rotation to the molecular frame and transformation 
to the LS basis states, 



E 



LS LS r , 0,0' o' 



^L'S'W^n'g 



X F] 



LSUl 



l]^ /sin0d0d0A^4,,j,i?-<.,o, 



1I* AT. J 



X {-i'L'S'n'^VL'sw'\di,\iLSVLLnsw) (B13) 
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where A). 



Dl .,i9,^,Q), and 



b labels the spherical basis components in the molecular 
frame. 

The terms involving J and J' can be expanded and the 
integration over the interatomic polar coordinates per- 
formed, to give 



where is the reduced matrix element of the dipole op- 
erator between the 2s and 2p atomic states. Only gerade 
{w = 0) ground states are coupled to the excited state, 
and because metastable states must satisfy (—1)'^"'"'" = 1 
due to their bosonic nature, only and states are 
coupled to the excited state. 



(2J' + 1)(2J+1) 



47r 



s\neded<t>Di,^^,D]:,D'^l,,^ 



2 J 



1 r-iJiJ' r^JiJ' 



(B14) 



dl between LS states 



The matrix element of dt, 
must be evaluated under proper symmetry considerations 
[2^. For the helium 2.s-2p(0+) transition this results in 



(L'(= l)S'n'i^n'sw\^ ^)\db\L{= 0)5f^L(= Q)nsw) 



V2 



(B15) 



After simplification the complete matrix element be- 
comes 



2 J' 4- 1 



^ ^n]'!-nj,o^mf,\,7ii'jd^ (B16) 



assuming that w' — 1 and w — 0. 
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